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PREFACE 


Turbulent  and  fluid -flow  phenomena  have  an  important  Influence  on 
atmoapheric  and  oceanic  circulation,  and  are  hence  expected  to  be  of 
considerable  concern  in  modeling  weather  situations,  and  estimating  the 
consequence  of  weather-modification  attempts.  For  reasons  of  simpli¬ 
fying  an  already  very  complicated  problem,  past  theoretical  turbulence 
work  has  been  largely  restricted  to  an  idealization  called  homogeneous 
turbulence.  In  this  report  an  effort  is  made  to  extend  one  kind  of 
turbulence  theory,  a  theory  for  "nearly  Gaussian"  processes,  to  the 
more  general  and  more  practical  class  of  problems  involving  nohhomogen- 
eous  flows. 

An  equilibrium  energy  spectrum  for  the  turbulence  in  such  prob¬ 
lems  is  found.  It  is  one  which  approximates  many  experiments  reason¬ 
ably  well.  A  hypothesis  for  one  term  in  the  theory  is  proposed,  and 
the  present  theory  is  used  to  correct  the  empirical  "mixing  length" 
theory  of  turbulence. 

Dr.  Meecham  performed  this  study  as  a  consultant  to  The  RAND 
Corporation.  He  and  Dr.  D.  -t.  Jeng  are  associated  with  The  University 
of  California,  Los  Angeles. 


ABSTRACT 


Earlier  work  with  the  nearly-normal  expansion  is  extended  to  in* 
elude  forced,  statistically  stationary  turbulence.  The  nature  of  the 
time-trans formation  of  the  Cameron-Martin-Wiener  representation  is  die- 
cussed.  For  a  stationary  process  we  obtain?* much  Information  using  the 
representation  at  one  time.  AFor  real  turbulence , we  find  E(k)  ~  k”  ; 
this  spectrum  is  discussed: and  shown  to  be  related  to  turbulence  inter* 
mittency.  Using  a  hypothesis  for  the  (small)  non-Gausslan  part  of  the 


velocity  --  a  hypothesis  suggested  by  earlier  work  --  a  corrected  mixing* 
length  theory  is  found  with  an  eddy  viscosity  rf^  vn_^L.u/i (a.  -  a9kI^V*) 
with  ai»a2  dimensionless,  positive  and  of  the  8ame"order;  k^  the 
wave  number  of  the  spectrum  maximum,  and  u '  and  £  the  rma  velocity 


fluctuation  and  the  turbulence  correlation  length  respectively. 
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I. _ INTRODUCTION 


In  recent  work  we  have  described  the  use  of  nearly-normsl  expan- 

(1  2) 

slons  (based  on  the  work  of  Cameron  and  Martin,  and  that  of  Wiener)'  * 
for  the  representation  of  turbulence.  The  expansions  were  used 

(3-5) 

to  deal  with  problems  of  decaying  turbulence,  n. ode  led  and  real. 

_2 

In  Ref.  3  the  correct  equilibrium  energy  spectrum  [E(k)  ~  k  ]  was  ob¬ 
tained  for  Burgers'  model  turbulence.  To  do  this,  a  hypothesis  was 
used  that  concerns  the  nature  of  the  non-Gausslan  term  In  the  velocity 
field.  The  same  result  can  be  obtained,  without  hypothesis,  for  a  sta¬ 
tionary  Burgers'  process. ^  It  Is  known  that  the  Cameron-Martln-Wlener 

/a:  q  \ 

expansion  is  not  uniform  in  time  for  nonforced  processes.  '  Thus, 
although  the  expansion  may  converge  rapidly  initially,  at  later  times 
the  convergence  becomes  less  rapid  even  for  nearly  normal  processes. 
Nevertheless,  It  has  been  found  that  the  convergence  Is  adequate  for 
the  treatment  of  real  wind-tunnel  turbulence,  for  up  to  six  correlation 
times  (or  120  mesh  lengths  downstream  from  the  wind-tunnel  grld).^ 

The  expansion  has  certain  virtues  not  possessed  by  some  other  represen¬ 
tations  for  turbulence.  These  advantages  have  been  emphasized  In  earlier 
work  and  will  not  be  recalled  here. 

In  most  practical  problems  involving  turbulence,  one  is  faced  with 
statistically  Inhomogeneous  turbulence,  which  Is  frequently  driven  and 
approximately  In  statistical  equilibrium.  The  driving  source  Is  very 
often  found  in  large-scale  effects  deriving  their  energy  from  the  mean 
flow.  For  example,  we  point  to  the  existence  of  turbulence  In  the 
lower  atmosphere  which  is  driven  by  wind  shear  in  that  region.  Another 
important  application  is  in  turbulent  wakes,  for  Instance  of  re-enter¬ 
ing  vehicles.  In  such  applications  we  find  the  turbulence  driven  by 
large  eddies  which  couple  into  the  mean  shear  flow  and  derive  their 
energy  from  that  flow. 

We  shall  here  replace  the  large-scale  energy  sources  by  an  equiva¬ 
lent  forcing  term.  We  suppose  that  the  process  has  existed  long  enough 
that  it  is  statistically  stationary  (in  a  statistically  inhomogeneous 
problem  the  forcing  term  may  vary  with  position).  It  Is  easy  to  gener¬ 
alize  the  random  process  expansion  for  such  problems  (see  Ref.  5,  Section 
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II).  Recently,  Saffman^^  has  used  this  kind  of  stationary  expansion 
to  deal  with  turbulent  diffusion,  and  has  obtained  sortie  very  promising 
results  in  this  way.  We  formulate  the  problem  for  stationary,  incom¬ 
pressible  real  turbulence  and  obtain  the  equilibrium  spectrum  in  that 
case.  The  result  is  k~  in  the  inertial  range.  Drydenv  '  has  observed 
approximately  such  spectral  behavior  for  turbulence  at  moderate  Reynolds' 
number.  This  spectrum  has  been  suggested  by  Townsend,  using  a  klnemat- 
ical  argument  (the  discussion  here  gives  the  results  using  the  dynamics). 

The  flow  realizations  leading  to  various  spectral  laws  will  be  dis- 

_2 

cussed.  In  particular,  as  is  known,  k  is  characteristic  of  flows 
with  near  discontinuities,  or  intermittences. 

A  hypothesis  is  proposed  concerning  the  form  of  the  non-Gaussian 
term  in  the  velocity-field  representation.  From  this  hypothesis  the 
inertia)  transfer  is  obtained  and  a  correction  to  mixing- length  theory 
is  proposed. 
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II .  DRIVEN  TURBULENCE 


We  develop  Che  equations  for  stationary  turbulence  in  an  lncom- 

(3-5) 

pressible  fluid,  relying  on  previous  work  for  some  of  the  details. 
Define  the  Fourier  transform  of  the  velocity  field 


t^Ck.t)  -  J  e1-*-  u1(r,t)dr 

(and  similarly  for  f).  The  Fourier  transform  of  the  incompressible 
Navier-Stokes  equation  is 

(ft  +  vk2)ut (IS.)  “  \  (2tt)"3  pij£^  J  ujQl  "  lc/)ujt(k/)dk/  +  fj^k.t) 

(2.1) 


with 


kiui 


0 


and 


Pijje(k)  "  kjgPiJ^  +  kjPU(k) 


and 


"  k’2ktkj  . 

The  given  forcing  term  £  will  be  assumed  to  be  statistically  homogen¬ 
eous  and  Gaussian  and  confined  to  large-scale  effects.  That  is,  it  is 
supposed  that  f  vanishes  for  k  larger  than  an  energy-containing  wave 
number  kQ.  The  first  two  random  processes  in  the  Cameron-Martin-Wiener 
representation  are 


4- 


tK^l)(r,ti) 


t"lj f  (— 1  ’—2  ’  W 


t"ll><-1,t1)  t^^’V  ‘  6ij  6(-l  ’  l2)6(tl  *  t2) 

(2.2) 


with  the  covariance  property 

<  tHll>(£l»tl)  •  6lj  6(^1  *  I2)6(tl  •  t2)  • 

The  presubscript  t  Indicates  a  transformation  (as  yet  unspecified)  In 
the  basic  whlte-nolse  process.  We  discuss  these  transformations  below. 

We  shall  suppose  later  that  the  turbulence  Is  statistically  sta¬ 
tionary  (It  Is  known  that  certain  necessary  conditions  must  be  Imposed 
for  statlonarlty);^12^  and  that  It  Is  (  t  least  locally)  statistically  ho¬ 
mogeneous  and  Isotropic.  We  suppose  eon  ’ergence  sufficiently  rapid 
that  the  random  field  can  be  adequately  represented  by  Just  the  first 
two  terms  In  the  nearly  normal  expansion  at  a  given  Instant.  Arguments 
have  been  presented  for  the  near  Gausslanlty  of  energy-range  turbulence. 

We  suppose  the  process  to  be  nearly  normal  In  the  Inertial  range  as 
well.  We  have 


ut(r,t)  a  u<l)  +  u<2) 


JJ  Ki(i)(-  *  t,t')  tH»l)(-/,c#)d-,dt# 
+  J/f J  K^0 ■  -l •  -2! 


i  (2) 


X  (r^r^;  ld— 2dtldt2 


(2.3) 
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^(r.t)  -  JJ  F^}d  -  r/;  t,t')  fcH^^(r /ft/)dr#  dt# 

+  JJJJf£b  -f  -  •  s*8  '-w 

X  1*— 2  *  ti»t2^d-ld-2dtldt2  <2**) 

We  suppose  that  our  representation  la  such  that  ve  can  use  the  manl- 
festly  statistically  homogeneous  form  involving  kernels  that  are  spa* 
tlal  functions  of  the  difference  variable.  The  transformation  of  the 
basis  in  time  represented  by  is  assumed  to  have  this  character* 

istic. 

Consider  the  behavior  of  the  Cameron-Mar tln-Wiener  (C-M-W)  representation 
(6-8} 

in  time.  It  is  known v  '  that  the  relation  of  the  random  process  to 
a  given  representation  changes  in  time,  in  many  nonlinear  dynamical 
problems ,  in  a  way  analogous  to  a  relative  rotation  of  a  vector  with 
respect  to  its  coordinate  system.  Thus  if  ve  consider  a  nearly  normal 
process,  u^,  at  some  instant  of  time,  ve  may  think  of  the  first  term 
in  the  representation  [i.e.,  the  first  term  of  (2.3)]  as  the  projection 
of  the  process  on  the  "Gaussian"  axis,  and  similarly  for  the  second 
term  on  the  lowest-order  "non -Gaussian"  axis.  Then  as  time  goes  on, 
although  the  process  may  remain  nearly  normal,  there  is  a  rotation  of 
the  process  relative  to  the  axes,  giving  different  components  at  later 
times  on  the  various  axes.  In  this  language  the  averages  (the  moments) 
correspond  to  vector -dot  products  of  the  various  quantities,  and  these 
products  are  Independent  of  the  coordinate  system  used  in  their  calcu¬ 
lation,  The  moments  are  correctly  given  by  any  representation  that  is 
complete,  and  the  one  used  here  is  complete,  so  ve  can  always  expect 
correct  results,  even  at  later  times  when  the  process  has  transformed. 

But  of  course  it  may  be  necessary  to  use  more  terms  in  the  representa¬ 
tion  as  time  goes  on.  There  is  a  "best"  representation  for  a  nearly 
normal  process  at  any  time,  which  yields  a  maximum  value  for  the  pro¬ 
jection  on  the  "Gaussian"  axis,  but  there  are  special  difficulties  for 
an  effort  directed  toward  finding  this  "best"  representation. 


This 
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klnd  of  transformation  of  the  representation  can  occur  even  for  sta¬ 
tionary  processes.  Thus,  although  at  a  given  time  the  stationary  pro¬ 
cess  may  be  representable  by  the  first  term  in  a  nearly  normal  expan¬ 
sion,  in  general  at  a  later  time  higher-order  terms  will  be  needed. 
Furthermore  for  the  specified  stationary  Gaussian  force,  £,  we  could 
at  any  Instant  use  one  term;  but  the  "best"  velocity-field  representa¬ 
tion  is  transforming,  and  since  we  must  use  the  same  representation 
for  the  two,  we  must  allow  the  possibility  of  transformation  even  for 
the  force.  This  reflects  in  the  requirement  that  the  kernel  for  the 
force  in  (2.4)  must  be  a  separate  function  of  the  times  t  and  t7,  and 
not  a  function  of  Just  the  difference.  At  a  given  time,  say  t  *  0, 
we  can  for  a  Gaussian  force  and  a  homogeneous  and  stationary  process, 
choose  a  "best"  representation.  The  force  then  has  only  the  first  term 
of  (2.4)  and  all  kernels  are  functions  of  the  space  difference  varia¬ 
bles.  Our  work  here  will  be  mainly  concerned  with  the  use  of  such  a 
fixed  representation. 

Returning  to  (2.1)  we  construct  the  usual  energy  equation.  Mul¬ 
tiply  by  the  velocity  at  a  second  wavenumber  (and  the  same  time),  sym¬ 
metrize  and  average  to  find, 

rfe  +  v<k*  +  k*>]  |  (2-0-3  P^Ck,) 

x  J  “  k#)Up(k#)uj(k2>)dk/ 

+  <f1(k1)Uj(k2)> 

+  same,  interchange  subscripts  l,j 
and  variables  k^,  k2 

(2.5) 


For  homogeneous  processes  we  know  that 

<ui (— ! )uj (— 2»  ~  V,  ki  +  k2  -  ° 


0,  otherwise 


(2.6) 
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where  V  is  the  volume  of  the  turbulence.  There  are  similar  relations 
for  other  moments. 

(14) 

Using  the  usual  definitions  we  have 

E(k)  -  T (2tt)2vJ  1  k2<Ulvk)u*(k)>  (2.7) 

and  for  isotropic  processes  (2.6)  becomes 

(|^  +  2vk2)  E(k)  -  T(k)  +  S(k)  (2.8) 

where 

T(k)  «  (2n)’5  k2Plcy0(k)  J  (u^k  -  k')u0(k (-k))dk/  +  c.c.  (2.9) 

S (k)  *  (f^ (k)u*(k))  +  c.c.  (2.10) 

(Actually  E  vanishes  for  stationary  processes.)  Here  c.c.  stands  for 
the  complex  conjugate  of  a  preceeding  term. 

From  the  properties  of  the  representation  (2.3)  we  know  that  to 
lowest  order  in  u  ,  the  transfer  can  be  written 

T(k)  -  y  (2n)'5k2PlaB(k)  J  <2u£l)(k  -  kOu^^k'^^-k) 

(2.U) 

+  u(1)(k  -  k/)u^1^(k,)u^2^(-k)>dk/  +  c.c. 

where  u^(k)  are  the  Fourier  transforms  of  u^(r),  time  dependence 

implicit.  To  write  the  transfer,  we  need  information  about  the  non- 

(2) 

Gaussian  term.  A  hypothesis  for  u  is  discussed  below.  T(k)  ■  0 
in  the  inertial  range,  i.e.,  for  k^  <  k  <  viscous  wave  number,  when 
the  process  is  in  equilibrium. 

For  a  stationary  process  T(k)  ■  0  for  all  k.  This  gives  us  a  con¬ 
dition  on  the  (dominant)  Gaussian  term  alone.  If  we  take  the  derivation 
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of  (2.9)  and  use  (2.1)  to  replace  u,  we  find  T  is  composed  of  fourth- 
order  moments  of  the  velocity  plus  third -order  moments  involving  the 
velocity  (alone  and  with  the  force).  Then  we  suppose  that  !u<2>|  «  |U(1>|, 
and  retain  in  the  fourth-order  moment  only  u^  terms  to  find  [correct 
to  terms  0(u  (2)2)-j 


T(k)  -  0  3? 


T(1)(k)  +  T(v)(k)  +  T(s)(k) 


(2.12) 


•  (v) 

where  T  stands  for  the  terms  proportional  to  v;  we  have 

T(1)(k)  -  Jk2  J  72  {2fk2(k"2  -  3k'2)  -  (k"2  -  k'2)2]$(1)(k,k',k") 


k  k  k 


+  2k"2 (k"2  -  k2)‘(1)(k',k",k)^dk' 


(2.13) 


where 


,n  E  (k")  [E  (k')  E  (k)1 

i<l)(k,k',k")  -  -^y-y  M-y - *5- 

16ttV2  Ik'2  k2  J 


(2.14) 


and 


Q  -  k4  +  k'4  +  k'4  -  2k2k'2  -  2k 'V2  -  2k,/2k2  (2.15) 

k  +  k'  +  k"  -  0  (2.16) 

•  ( 1  \  •  (v) 

The  time  derivative,  T  +  T'  ,  of  the  transfer  is  exactly  that 
obtained  by  Proudman  and  Reid^15^  using  the  zero-fourth-cumulant  hy¬ 
pothesis.  (We  employ  a  notation  similar  to  that  of  Ref.  15  here.  The 
viscous  term,  T^V\  will  not  be  used  here.) 

In  (2.13),  Ej  is  the  "Gaussian"  part  of  the  energy.  It  is  obtained 


from  (2.7), 
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Ex(k)  -  f2rT2v]  1k2<u^l)(k)u^l)(-k». 


(2.17) 


From  (2.2)  and  (2.3)  using  the  transform  we  have  at  the  initial 
instant  t  =  0, 


E1  (k,0)  =  2(2rr)"2k2  J  |K(1)(k;  0ft')|2dt'  (2.18) 


with  defined  implicitly  by 

K^^k;  0,t')  =  PiJ(k)K(1)(k;  0,t')  (2.19) 

and  with  [see  (2.3)") 

K^}(k;  0,t')  -  j  elk*—  K^j}(r;  0,t')dr  .  (2.20) 

To  repeat,  it  is  understood  that  the  representation  will  change  at  a 
later  time;  we  use  the  representation  at  t  ■  0. 

The  source  contribution  is 


T(S)  =  y  (2-)"5  k2Pi(y0(k)  J  dk'  <2u^(k  -  k/)fg(k')ui(-k) 


+  u  (k  -  k;)u  (k,)f.(-k))  +  c.c. 
a  —  —  9  -  i  - 


(2.21) 


* (S )  ( 2 ) 

To  lowest  order  each  term  for  T  initially  is  of  order  u  Just  as 
for  the  transfer  (2,11). 

There  are  important  differences  between  the  work  in  Ref.  15  and 
the  present  discussion.  The  process  here  is  statistically  stationary 
(and  of  course  driven).  The  equations  can  be  used  only  at  or  near  a 
nearly  normal  statistical  equilibrium.  In  order  to  follow  the  approach 
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to  equilibrium,  we  should  have  to  integrate  the  equations  for  the  ker¬ 
nels  in  the  nearly  normal  expansion  in  the  usual  way.^*^  If  the  ini¬ 
tial  process  is  near  enough  to  its  equilibrium  form,  this  can  be  done 
with  a  rapidly  converging  series.  However,  if  the  initial  process  is 
far  from  equilibrium,  the  intermediate  stages  may  be  highly  non-Gaus- 
slan,  and  the  series  poorly  convergent.  Then  the  Gaussian  part  of  the 
derivative  of  the  transfer,  of  (2.13),  will  not  be  a  good  approx¬ 

imation  and  the  expansion  falls.  This  has  occurred  in  some  previous 
computations  making  use  of  the  zero-fourth-cumulant  approximation, (^,18) 

where  (nonphysical)  negative  energy  spectra  have  been  found. 

•(v) 

To  determine  T  ,  we  note  that  we  need  the  Fourier  transforms  of 

the  triple  velocity  correlation.  A  knowledge  of  the  non-Gaussian 

(2) 

term  u  in  the  viscous  range  is  needed  for  their  evaluation.  We 
shall  examine  only  the  energy  and  inertial  ranges  in  this  paper. 

The  time  derivatives  of  all  moments  vanish  for  stationary  processes. 
One  may  wonder  why  we  focus  attention  on  T(k)  [see  (2.12)].  In  the  in¬ 
ertial  range  the  derivatives  of  even  moments  are  small  for  nearly  normal 

(2) 

processes  since  those  derivatives  involve  u  .  Tne  derivative,  T,  In¬ 
volving  the  first  important  odd  moment  is  chosen  because  of  its  central, 
physical  role. 
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III.  INERTIAL  RANGE  EQUILIBRIUM 

At  equilibrium  we  have  in  the  inertial  range  [see  (2.13)],  assum¬ 
ing  we  can  separate  the  viscous,  inertial,  and  energy  ranges  in  the 
usual  way, 

0  «  T(1)(k)  -  J  k2  J  -2“^2  {2[k2<k"2  -  3k/2>  "  <k"2  -  k,2)2]*(k,k/,k/') 

+  2k/,2(k//2  -  k2)$(k',k\k)}dk'  (3.1) 


with 


*(k,k',k") 


E^k")  p1(k/)  E^k) 

i  c.  2«_  // 2  .  /2  .2 

16tt  k  L  k  k 


(3.2) 


for  kn  <  k  <  k  ,  and  k  is  that  wavenumber  at  which  viscous  effects  are 
0  v’  v 

no  longer  negligible. 

We  continue  by  solving  for  the  spectrum  E^,  which  will  cause  the 

derivative  of  the  transfer  to  vanish  in  the  inertial  range.  Note  first 

2 

the  interesting  result  that  the  "equipartition"  form  E  ~  k  is  an 
(exactly  Gaussian)  solution.  Look  for  an  algebraic  form,  E^  ■  (c/2rr2)k”n. 
Substitute  in  (3.1)  -  (3.2),  make  changes  of  variable,  and  find 


0  - 


(3.3) 


with 


I 


j 


-n , 
r  dr 


pCj  (1  -  u2)P1(r,n) 

J  -1  a\ 
(l  -  2ru  +  r2) 


du, , 


J  "  1,2,3 


(3.4) 


with  quantities  given  In  Table  1 
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Table  1 

QUANTITIES  IN  (3.4) 


We  look  for  the  root  of  (3.3);  the  right  side  diverges  (with  op- 

■f  .  • 

poslte  sign)  for  n  -  1  ,  3  .  A  quantity  proportion  to  T  —  the  right 
side  of  (3.3)  --  is  plotted  In  Fig.  1.  The  calculation  Is  complicated 
by  cancelling  singularities  in  (3.4).  The  root  is  clearly  very  near 
to  n  ■  2,  If  not  exactly  that  value.  We  have  not  calculated  the  dashed 
portion  of  the  curve  In  Fig.  1;  we  merely  sketch  those  divergent  por¬ 
tions  . 

m2 

It  Is  known  (see  Ref.  3)  that,  for  a  one-dimensional  process,  k~ 
is  the  spectrum  which  is  characteristic  of  near -discontinuities  in  the 
flow.  It  is  easy  to  show  that  this  is  also  the  case  for  some  three- 
dimensional  flows  (see  Appendix).  Real  turbulent  flows  certainly  ex¬ 
hibit  near -discontinuities  and  might  be  supposed  to  have  an  energy  spec- 
trum  range  with  a  k  behavior.  Indeed  early  measurements  showed  ap¬ 
proximately  this  behavior,  as  can  be  seen  in  Fig.  2. 

•2 

From  (3.1)  we  see  we  cannot  determine  the  coefficient  of  k  . 

Guided  by  earlier  experimental  works  (Fig.  2  and  Ref.  14,  Fig.  7.8) 
for  the  energy  range  we  write 

E(k)  -  C  U,2k~2 

I 


(3.5) 
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U 


Fig.  2  —  Energy  spectra  from  Dryden,  Ref.  13. 


Comparison  of  National  Bureau  of  Standards  and  National  Physi¬ 
cal  Laboratory  measurements  of  the  spectrum  of  turbulence, 
plotted  nondimensionally. 

At  left,  NBS  values  40  (•)  and  160  (+)  Inches  behind  1-lr.ch 
mesh  screen  at  U  =  40  ft/sec. 

At  right,  NPL  values  82  inches  behind  3-inch  mesh  screen  at 
15  (•),  20  (x),  25  (+),  30  (L) ,  and  35  ft/sec. 

The  reference  curve  in  each  case  is  the  curve 


UF (n)  B 


,22.2 

4tt  n  ly 


1  + 


U 


00 

U  is  the  mean  speed,  L  is  the  integral  J  R  dx,  R  is  the 

o 

correlation  between  the  fluctuations  at  two  points  separated 
by  the  distance  x  in  the  direction  of  flow,  n  is  the  frequency 
and  E(n)  is  the  fraction  of  the  total  energy  of  the  turbulence 
arising  from  frequencies  between  n  and  n  +  dn. 
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e  =  J  T(k)dk  (3.6) 

Energy 

Range 

with 

€  =  |  F.(k)dk  *  (3.7) 

The  Kolmogoroff  spectrum  ^k  may  depend  upon  higher-order 

terms  In  the  C-H-W  series. 
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IV .  MIXING-LENGTH  CORRECTION 

To  find  T(k)  for  the  statistically  stationary  process  given  in 

(2) 

(2.11)  we  see  from  (2.3)  that  we  shall  need  a  knowledge  of  K  given 

in  terms  of  the  representation  at  t  ■  0,  i.e.,  in  terms  of  _H^. 

(5)  a 

In  calculating  for  turbulence  experiments,  good  results  have  been 
obtained  using  a  hypothesis  for  this  kernel  (stated  in  terms  of  the 
Fourier  transform) 


Kijk%,!52*  "  2  C0  Piy6(-1  + -2)KYj>(-l)K6k>(-2)#  (4,1) 

(2) 

Incidentally  this  is  the  form  of  K  for  short  times  if  the  process 
is  initially  exactly  Gaussian  [see  Ref.  5,  Eq.  (3.14)].  Here  tQ  is  a 
characteristic  tine  of  order  the  correlation  time;  in  Ref.  5  t0  (there 
called  Oq)  was  found  to  be  0.4  for  some  turbulence  experiments.  We 
substitute  (4.1)  in  (2,11)  and  use  averages  of  products  of  processes 
(see  Ref.  5)  and  (2.18)  to  find 


T(k) 


k/2(k.k')  +  2k2(k*k/)] 


y  EM  EXkp.dk' 


k' 


2 


r  Q(k,k/,k//)  r 

k2k/2k//2 


-  k2(k'*k") 


2(k*k')(k*k")] 


X 


E(k')  E(k") 
k/2  k"2 


(4.2) 


We  use  the  definitions  of  Section  II  and  assume  that  the  energy  spec> 
trum  is  primarily  given  by  its  Gaussian  part  [here  Just  E(k)]  at  least 
for  moderate  wavenumbers. 

It  is  interesting  that  this  transfer  is,  except  for  a  coeffici- 
•  (1 )  * 

ent,  equal  to  T  ,  the  part  of  T  using  lowest  order  (Gaussian) 


-17- 


contributions  [see  (2.12)  through  (2.15)].  In  the  Inertial  range, 

-2 

E  ~  k  gives  zero  transfer  using  (4.1)  --  as,  in  the  previous  section 
It  gave  zero  T.  Thus  (4.1)  is  one  solution  for  giving  T  *  0  in 

the  inertial  range.  The  same  correspondence  held  for  Burger's  model 
(see  Ref 8.  3  and  6).  We  are  Interested  in  this  section  in  E  and  T 
for  small  k  in  the  energy  range. 

Consider  a  version  of  the  mixing-length  theory.  One  can  define 
an  eddy  viscosity  to  account  for  the  nonlinear  transfer  [refer  to  (2.8)], 
treating  the  effect  of  that  transfer  on  the  large  eddies, 

2vD(k)k2E(k)  =  -T(k)  (4.3) 


or 


vD(k)  =  -T(k)/[2k2E(k)]  .  (4.4) 

It  is  sometimes  convenient  to  use  (4.3)  with  vD  approximated,  in 
place  of  the  full,  nonlinear  problem.  We  expand  the  expression  for 
vD  about  k  =  0  using  the  functions  above. 

To  lowest  order,  the  angular  integrations  in  (4.2)  for  T(k)  can 
be  performed.  The  result  of  these  manipulations  is 

VD(0)  =  I?  t0  {2ET  “  E" (0)  /2  ^  dk/}  * 

v  K 

We  suppose  that  the  curvature  of  E(k)  doesn't  vanish  at  the  origin  for 
the  nonhomogeneous  turbulence  being  approximated  here.  Also  we  define 
the  turbulent  energy  per  unit  mass, 


E 


T 


E(k')dk' 

0 


(4.6) 
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The  second  term  in  (4.5)  tends  to  emphasize  the  larger-scale  part  of 
the  process,  reducing  the  eddy  viscosity.  The  eddy  viscosity  given 
by  (4.5)  is  clearly  of  the  order  of  the  usually  assumed  value,  for  we 
have 


E 


T 


u 


/2 


with  u'  the  fluctuation  velocity,  so  [from  (4.5)] 


VD  ~  C0  ET 


u  't 


(4.7) 


with  j l  the  turbulence  scale.  In  the  theory  under  discussion,  it  is  not 
evident  from  the  start  that  a  form  like  (4.5)  would  be  obtained  at  k  ■  0. 
The  result  is  felt  to  be  encouraging  for  the  use  of  the  nearly  normal 
expansion. 

We  now  continue  to  the  next  order  term  in  k.  To  simplify  the  dis¬ 
cussion  we  use  an  estimate  for  E(k).  This  restriction  isn't  necessary 
but  is  convenient.  Set 


E(k)  -  \  E"(0)k2 


6__  .2/3  -5/3 

n  ET  E  * 


k  <  k. 


k  >  k. 


(4.8) 


where  k^  is  the  position  of  the  maximum  of  E,  and  E  is  assumed  to  be 
continuous  there.  The  larger  wavenumber  behavior  is  of  course  the  Kol- 
mogoroff  spectrum.  The  result,  to  lower  order,  for  (4,4)  is 


vD(k) 


tQ  ET  [o. 123  +  0.139 


s  vD(0)  +  \  v^(0)k2 


(4.9) 
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Use  of  the  spectrum  (4.8)  permits  us  to  calculate  all  Integrals  in¬ 
volved.  We  can  plausibly  replace  the  nonlinear  terms  in  the  Navier- 

2 

Stokes  equation  by  v^(k)k  u.  In  physical  space,  the  modified  Navier- 
Stokes  equation  is 


0  bt  =  {VD(0)  ’  \  '^0)*2  +  4  '^u  *  (4*10) 

where  in  turbulent  flows  v  is  negligible.  The  coefficients  are  found 
from  (4.9).  (It  is  interesting  that  the  two  coefficients  in  (4.9)  are 
of  the  same  order.)  The  relation  provides  a  one-parameter  form  for  the 
two  coefficients  (remembering  that  k^  is  the  position  of  the  spectral 
maximum).  An  experimental  check  would  be  helpful. 


-20- 


V.  CONCLUSIONS 


We  have  considered  the  use  of  the  Cameron-Martin-Wiener  represen¬ 
tation  for  statistically  stationary  (and  inhomogeneous)  turbulence. 
Difficulties  with  the  time  transformation  of  the  white-noise  process 
are  avoided  by  using  the  representation  only  at  one  instant.  By  deal¬ 
ing  with  the  time  derivative  of  the  energy-transfer  term,  we  obtain  a 
formulation  similar  to  the  zero-fourth-cumulant  approximation  --  though 
differing  from  it  in  important  ways.  From  this  formulation,  a  k*  en¬ 
ergy  spectrum  is  obtained  and  discussed.  Using  a  hypothesis,  supported 
by  previous  computations,  a  useful  correction  to  mixing-length  theory 
is  found. 
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Appendix 

FLOW  REALIZATIONS  WITH  E  ~  k 


Some  velocity  flows  with  near-discontinuities  have  spectra  with 
large  wavenumber  behavior  like  k  Townsend  has  shown  this  spec¬ 
trum  results  for  flows  with  random  shear  discontinuities,  i.e.,  vortex 
sheets.  Consider  a  flow  made  up  of  spheres  rotating  with  angular  ve¬ 
locity  a).  We  examine  first  the  Fourier  transform  of  one  such  sphere: 

u(k)  =  |JJ  e1-*-  u(r)dr  (A.l) 


with  (in  spherical  coordinates) 

A 

u(r)  =  0o)r  sin  0,  r  <  a 
“0,  r  >  a 


(A.  2) 


where  $  is  the 

a 

e.g.,  uy(ki). 
that  for  ka  » 


unit  azimuthal  vector.  Let  us  consider  a  typical  value, 
Using  the  method  of  stationary  phase  it  is  easy  to  show 
1, 


uy(ki) 


(A. 3) 


Suppose  we  have  a  random  set  of  N  such  spheres  in  a  finite  volume 

(14) 

V.  The  energy  spectrum  is  found  (using  volume  averaging) 


E(k) 


2  2 

k  Nu  (k) 
V 


(A.4) 


From  (A. 3)  we  see  a  typical  value 


E (k)  ~ 


M  2  * 

Nu;  a 


V 


(A.  5) 
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It  is  easy  to  show  many  other  discontinuous  flows  with  this  spectral 

behavior.  It  is  known  that  the  corresponding  behavior  holds  for  nearly 

(3) 

discontinuous,  one-dimensional  time  series. 
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